Finite temperature theory of the trapped two dimensional Bose gas 
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We present a Hartree-Fock-Bogoliubov (HFB) theoretical treatment of the two-dimensional 
trapped Bose gas and indicate how semiclassical approximations to this and other formalisms have 
lead to confusion. We numerically obtain results for the quantum mechanical HFB theory within 
the Popov approximation and show that the presence of the trap stabilizes the condensate against 
long wavelength fluctuations. These results are used to show where phase fluctuations lead to the 
formation of a quasicondensate. 
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The question of whether a weakly interacting Bose 
gas can undergo Bose-Einstein condensation (BEC) when 
confined to an effectively two dimensional (2D) geometry, 
has gained significant topical interest with the advent of 
recent experiments pj . It is well know that a 2D homoge- 
neous ideal gas of bosons does not undergo BEC at finite 
temperature Q- Indeed it has been rigorously shown 
that even when interactions between the particles are in- 
cluded there is no BEC phase transition at finite tem- 
perature, although a superfluid phase transition in the 
form of a Kosterlitz-Thouless (KT) transition |j| can be 
shown to take place. The KT transition occurs because 
of the enhanced importance of fluctuations in the 2D sys- 
tem yielding a quasi-condensate state where the length 
scale for phase coherence is small compared to the sys- 
tem size. The observation of such a phase transition has 
been reported recently p| for a 2D gas of dilute hydrogen 
on a liquid helium surface. For the non-interacting gas, 
Bagnato and Kleppner || showed that an ideal gas in 
a 2D harmonic trap does exhibit a BEC phase. Petrov 
et al. [l7| included interactions and showed within the 
Thomas-Fermi approximation, that well below T c there 
exists a true condensate while at higher temperatures a 
quasicondensate forms. We are interested in the phase di- 
agram at temperatures below and within this fluctuation 
regime to investigate both the true BEC state and the 
onset of phase fluctuations that lead to the destruction 
of the BEC. 

In three dimensions (3D), the finite tempera- 
ture Hartree-Fock-Bogoliubov (HFB) treatment of the 
trapped Bose gas has proved very successful Q, H, Hi- 
lt is therefore natural to develop this approach for the 
2D gas. Remarkably the quantum mechanical HFB for- 
malism has, to the best of our knowledge, never been 
implemented in 2D, the development of which is the cen- 
tral result of this article. First let us consider however, 
a simplification to the fully quantum mechanical imple- 
mentation of HFB that has been used to conclude pre- 
viously that there is no BEC in an interacting 2D Bose 
gas. 

In 2D, with a cylindrically symmetric trapping po- 
tential, the grand-canonical many-body Hamiltonian is 



given by 
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where we have assumed a spatially dependent coupling 
parameter, introduced by the many-body T-matrix for 
the trapped case as proposed in [l(J, [HI • 

In a manner identical to the development of the the- 
ory in 3D ^3] we assume that the many-body Bose 
field operator can be decomposed into a mean, conden- 
sate part and a fluctuating field operator part such that 
ip = (tp) + if) = 4> + tp. The above Hamiltonian can then 
be diagonalised provided the condensate order parame- 
ter, <fr, obeys the generalised Gross-Pitaevskii equation 
(GPE) 
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and the quasiparticle excitations obey the Bogoliubov-de 
Gennes (BdG) equations 
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where C = 
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is the condensate density and n the noncondensate den- 
sity, which is evaluated by populating the quasiparticle 
levels according to the usual Bose distribution. In the 
derivation of these equations we have used the canonical 
transformation ijj — • \uj&j — v*a\] where the &i satisfy 
the usual Bose commutation relations and we have taken 
the so called Popov approximation, neglecting anomalous 
pair averages of the fluctuating field operator. 

Together, the GPE and BdG equations form a closed 
set and can be solved numerically using techniques anal- 
ogous to those developed in 3D. This we will describe 
shortly, but first let us attempt to obtain a semiclassical 
solution of these equations. This is the approach taken 
by previous authors [l3| . 
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Let us assume that the temperature is large compared 
to the energy level spacing in the trap. We can now re- 
place the quasiparticle amplitudes such that Ui w ue l6 
and Vi rs ve lS where the common phase, (9, defines a 
quasiparticle momentum, p = HV9. Neglecting spatial 
derivatives of u, v, and p (local density approximation) 
and making a continuum approximation for the Bose dis- 
tribution function, the quasiparticle excitation spectrum 
is found to be given by 



E sc (p,r) = V - {gn c 
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condensate density is still calculated via the GPE. The 
density of excited state particles can now be integrated 
in closed form to give 
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equations form a closed set and it is this set of equa- 
tions that previous authors have attempted to solve self- 
consistently. This is not possible however, as the semi- 
classical approximation can only be used consistently at 
low energy if combined with a Thomas-Fermi approxi- 
mation for the condensate [T^j . There are therefore no 
solutions to these equations. 

It is trivial to show that in the Thomas-Fermi limit 
t = s and the expression for the semiclassical thermal 
density is undefined. Indeed, even for low particle num- 
bers, it can be shown that the arguments of the square 
roots always contain values at some spatial point that are 
negative (and approach zero from below in the Thomas- 
Fermi limit) and hence this expression is never well de- 
fined. The origin of this problem lies in the expression for 
the semiclassical excitation spectrum. At low energies, or 
equivalently long wavelengths, Eq. (4) yields imaginary 
energies. This has been used to conclude that in 2D BEC 
cannot take place since the condensate is destabilised by 
long wavelength fluctuations |l-3|. This is simply incor- 
rect. What one is seeing is a failure of the semiclassical 
approximation. If the discrete nature of the excitation 
spectrum for the finite size trap is not retained for the 
low energy excitations, and this is the case within the 
semiclassical approximation, there comes a point where 
the magnitude of the chemical potential exceeds the mag- 
nitude of the effective potential and the argument in Eq. 
(4) becomes negative for low p. Just because the semi- 
classical treatment of the HFB formalism fails does not 
mean that the quantum mechanical theory will also. In 
this case the long wavelength oscillations, corresponding 
to the p = limit in the semiclassical case, are precluded 
by the finite size of the trap. Therefore, to determine 
whether BEC can take place in 2D, the full, numerically 
expensive, discrete calculation must be undertaken. 

As an aside, it is possible to obtain a well defined semi- 
classical theory of the trapped 2D gas if one makes the 



further Hartree-Fock approximation [l5| . This consists of 
setting the v(r, p) terms in the semiclassical HFB treat- 
ment to zero everywhere. The Hartree-Fock excitation 
spectrum is now single particle-like and is just given by 
A sc rather than by the phonon-like Bogoliubov spectrum 
of Eq. (4). There is therefore no problem with the in- 
frared divergence seen previously and one can integrate 
to obtain the semiclassical Hartree-Fock thermal density, 
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which can be solved self-consistently together with the 
GPE to obtain a Bose condensed solution. Interestingly, 
if one omits the ground state from the calculation, ef- 
fectively demanding that condensation does not occur, 
then it is still possible to obtain a self-consistent solution 
for the thermal density at all temperatures. Therefore, 
at the level of the semiclassical Hartree-Fock approxi- 
mation, one might conclude that there is no BEC in the 
thermodynamic limit for a 2D gas of trapped, interacting 
bosons. 

By removing the lowest lying state however, one has 
truncated the available state space and so one may ex- 
pect that the thermodynamically stable state is the con- 
densed phase. This is confirmed by calculating the free 
energies for the two states. Above the BEC transition 
temperature the free energy of the condensed and un- 
condensed solutions are virtually indistinguishable, but 
below the critical temperature the free energy for the 
condensed solution is the lower, confirming that within 
the Hartree-Fock approximation, BEC is the thermody- 
namically favoured state |13|. 

We now proceed to solve the quantum mechanical, 
discrete, HFB equations [Eq. (2) and Eq. (3)] self- 
consistently. The method of solution is described in de- 
tail elsewhere |l2|. but is outlined here for completeness. 

First, the GPE equation is solved using an expansion 
in some appropriate basis set with the condensate and 
noncondensate densities set to zero. The solution for the 
condensate wavefunction is then used to construct n c and 
the process iterated to find a converged solution. One 
now needs to calculate the noncondensate density. To 
do this one decouples the BdG equations by making the 
transformation to the auxiliary functions ip^ = itj ± 
Vi. One can thus obtain equations for tp^ and tp\ ^ 
separately. These are solved using a further basis set 
expansion, for which we use the basis of excited states of 
the GPE to ensure orthogonality with the condensate and 
to simplify the construction of the matrix elements. The 
noncondensate density is then constructed by populating 
the quasiparticle states. This value of n is inserted into 
the GPE and the whole process repeated iteratively until 
convergence. 

In contrast to the semiclassical HFB treatment we 
have no difficulty in finding self-consistent solutions to 
the quantum, discrete HFB equations and typical results 
are presented below. In the numerical solution we non- 
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FIG. 1: Quasiparticle excitation frequencies as a function of 
temperature for N = 2000 atoms and the full spatially de- 
pendent coupling parameter. 
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FIG. 2: Fraction of atoms in the condensate as a function 
of temperature. The dotted line corresponds to the non- 
interacting gas for comparison. 



dimensionalise our equations using the natural harmonic 
oscillator units. In this case we take our Rydberg of en- 
ergy to be Kujq/2 and in these units the non-dimensional 
interaction parameter, g(r), takes values between 0.09 
and 0.1. This parameter, in 3D, depends only upon the 
s-wave scattering length, a, however, in 2D, g also de- 
pends upon the strength of the confining potential in the 
third dimension. The 2D gas thus represents a system 
where the interparticle interaction strength can be tuned 
by modulating the confining potential in the z-direction 
in a manner analogous to the use of the Feshbach reso- 
nance in 3D [l^. It has even been suggested that the 
sign of the interaction strength (and hence whether the 
interactions are attractive or repulsive) can be changed 
by tuning the confining potential in the third direction 

In Fig. (1) we display the low lying excitation spectrum 
as a function of temperature for the trapped gas of 2000 
87 Rb atoms. Shown are the lowest lying m = 1, m = 
2 and to = quasiparticle, or collective, modes up to 
the critical temperature. We point out that the results 
shown are those of the HFB-Popov equations. In three 
dimensions, in the presence of significant direct driving 
of the thermal cloud, the theory fails to predict shifts in 
the frequencies above about 0.6 T c These are related 
to the dynamics of the thermal cloud not included in this 
theory. The excitation spectrum shown should be valid 
in experiments where the condensate is driven directly. 

In the three dimensional asymmetric case this approx- 
imation is valid only below 0.6 T c . However, we include 
the higher temperatures for the symmetrical case where 
the theory in 3D does not fail. The m = 1 mode is the 
Kohn mode [19j and clearly satisfies the generalised Kohn 
theorem to within our resolution, apart from a slight de- 
viation near the critical temperature. The quadrupole 
(m = 2) and breathing (to = 0) modes are well defined 
and nowhere look like going soft (approaching zero fre- 
quency as an indication of an instability). We therefore 
conclude there are well defined solutions to the quantum 
mechanical HFB equations, in a manner exactly the same 
as in 3D, for all temperatures below the critical temper- 



ature. 

It is important to note that, at low temperatures, the 
frequency of the lowest lying to = mode, or breathing 
mode, is at precisely 2tt>o, independent of the interaction 
strength. This result was predicted some time ago by 
Pitaevskii and Rosch pfij. purely from symmetry argu- 
ments. Indeed, they show rigorously that the existence 
of 2u) oscillations is ensured by the underlying SO(2,l) 
symmetry of the full quantum theory Hamiltonian for the 
interacting harmonically confined 2D gas with a contact 
interaction. It is therefore noteworthy (and critical) that 
our calculations confirm this result numerically at low 
temperatures. At higher temperatures the result is still 
valid, however the excitations are those for a condensate 
in an effective potential which is modified by the addi- 
tion of the potential 2gn from the static thermal cloud. 
The condensate effectively sees a weaker harmonic po- 
tential and hence the to = mode has a slightly lower 
frequency. Above the critical temperature, the excitation 
frequencies, of course, go over to those of the thermal gas. 
Similarly, at high temperatures (near the critical temper- 
ature) the to = 1 mode no longer satisfies the Kohn the- 
orem precisely. This is again due to the presence of the 
potential from the thermal cloud. The effective potential 
in which the condensate oscillates is no longer harmonic 
and the Kohn theorem broken. As discussed previously 
[l2|. if the full dynamics of the thermal cloud were in- 
cluded, then the Kohn theorem would be satisfied at all 
temperatures. 

In Fig. |2 we show the condensate fraction as a func- 
tion of temperature. It is clear that at temperatures be- 
low the critical temperature we obtain a macroscopic oc- 
cupation of the ground state, which implies BEC. This 
however is not sufficient. For true condensation we re- 
quire a well defined phase over the entire condensate. In 
3D this is true for all temperatures below T c except for 
a small region near the critical temperature often ref- 
ered to as the Ginzburg region. In this region phase 
fluctuations prohibit the formation of a true conden- 
sate. In a uniform 2D gas this region extends all the 
way to T — and is what prevents the formation of a 
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FIG. 3: Single particle correlation function for the 2D Bose 
function of position at 0.05, 0.1, 0.35, 0.75 and 0.925 
T c from right to left, showing decreasing correlation length as 
a function of temperature. Lengths are scaled with the size 
of the condensate, r con , at each temperature. 

BEC. In Fig. we plot the off-diagonal correlation func- 
tion, g^'(0,r) |2l| showing that only at low tempera- 
tures is there a coherent condensate with a correlated 
phase spanning the condensate. As the temperature is 
increased the coherence begins to decay on a length scale 
less than the dimensions of the gas. In analogy with 
Petrov et al. 17], if we expand the field operator as 
ijj = 0o+<5"0 = Vn~e l< ^ then we can express the nonconden- 
sate density asii = (SffiSijj) — + i[5h, <j>]/2 + n </> 2 ). 
At low temperatures density fluctuations are suppressed 
and this yields n/no ~ (4> 2 )- Using the Thomas-Fermi 
approximation Petrov et al. conclude that phase fluc- 
tuations lead to the formation of a quasicondensate at 
temperatures of around T = T c /2 for our parameters. 
In our case we find that the condensate persists to ap- 
proximately this temperature, but above this the phase 



becomes ill defined and the system is best described as 
a quasicondensate. This is consistent with our numerical 
calculation of the correlation function. A treatment of 
the 2D gas which could describe the quasicondensate has 
recently been formulated |22J, but has yet to be imple- 
mented. 

In conclusion, we have shown that the presence of the 
trap stabilizes the condensate against long wavelength 
fluctuations. This is true not only for density fluctua- 
tions, but for phase fluctuations as well, which are in- 
cluded in our formalism via the contribution to the non- 
condensate density from low-energy quasiparticles. Our 
work is consistent with Petrov et al. 17]. A 2D trapped 
dilute gas of weakly interacting bosons therefore does un- 
dergo BEC, forming a pure condensate at temperatures 
below a transition region near T c . Although this conclu- 
sion has been reached by, among others, Petrov et al. 
and Bagnato and Kleppner the converse conclusion 
has also appeared in the literature 0]. The prospect of 
performing the HFB calculation has been proposed as a 
means of clarifying the issue (in addition to the above 
references, see for example, Bayindir and Tanatar |23|L 
but to this point no one had done so. We have now per- 
formed this calculation and unambiguously shown that 
BEC does occur for the 2D trapped interacting gas when 
the discrete nature of the energy spectrum is taken into 
consideration. 
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